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Abstract 

Two complementary methods to describe the collective motion, RPA and Wigner 
function moments method, are compared on an example of a simple model - har- 
monic oscillator with quadrupole-quadrupole residual interaction. It is shown that 
they give identical formulae for eigenfrequencies and transition probabilities of all 
collective excitations of the model including the scissors mode, which here is the 
subject of our special attention. The exact relation between the variables of the two 
methods and the respective dynamical equations is established. The normalization 
factor of the "synthetic" scissors state and its overlap with physical states are cal- 
culated analytically. The orthogonality of the spurious state to all physical states 
is proved rigorously. 



1 



1 Introduction 



The full analysis of the scissors mode in the framework of a solvable model (harmonic 
oscillator with quadrupole-quadrupole residual interaction (HO+QQ)) was given in [1]. 
Several points in the understanding of the nature of this mode were clarified: for example, 
its coexistence with the isovector giant quadrupole resonance (IVGQR), the decisive role 
of the Fermi surface deformation, and several things more. 

The Wigner Function Moments (WFM) method was applied to derive analytical ex- 
pressions for currents of both coexisting modes (for the first time), their excitation en- 
ergies, magnetic and electric transition probabilities. Our formulae for energies turned 
out to be identical with those derived by Hamamoto and Nazarewicz [2] in the frame- 
work of the RPA. This fact generated the natural motivation for this work: to investigate 
the relation between formulas for transition probabilities derived by two methods. More 
generally we will perform a systematic comparison of the two approaches. The HO+QQ 
model is a very convenient ground for this kind of research, because all results can be 
obtained analytically. There is no need to describe the merits of the RPA - they are very 
well known [3]. It is necessary, however, to say a few words about the WFM. Its idea is 
based on the virial theorems of Chandrasekhar and Lebovitz [4]. Instead of writing the 
equations of motion for microscopic amplitudes of particle-hole excitations (RPA), one 
writes the dynamical equations for various multipole phase space moments of a nucleus. 
This allows one to achieve a better physical interpretation of the studied phenomenon 
without going into its detailed microscopic structure. The WFM method was successfully 
applied to the study of isoscalar and isovector giant multipole resonances and low-lying 
collective modes of rotating and nonrotating nuclei with various realistic forces [5]. The 
results of WFM were always very close to similar results obtained with the help of RPA. 
In principle, this should be expected, because the basis of both the methods is the same: 
Time Dependent Hartree-Fock (TDHF) theory with its small amplitude approximation. 
On the other hand, they are equivalent only in special cases. The detailed analysis of the 
interplay of the two methods turns out to be useful also from a "practical" point of view: 
firstly, it allows one to obtain additional insight into the nature of the scissors mode; 
secondly, we find new exact mathematical results for the considered model. 

The paper is organized as follows. In Section 2 we remind the principal points of the 
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WFM formalism and give the summary of the key results of [1] obtained by applying this 
method to the HO+QQ model. The same model is considered in Section 3 in the frame 
of RPA: the formulae for eigenfrequencies, electric and magnetic transition probabilities 
of the scissors mode are derived, the "synthetic" scissors and spurious state are analyzed, 
the RPA equations of motion for transition matrix elements are compared with the WFM 
equations of motion for irreducible tensors. The exact relation between the RPA and 
WFM variables is established in Section 4. The mutual interplay of the two methods is 
discussed in the conclusion. The various mathematical details are given in Appendices A 
and B. 



2 The WFM method 



The basis of the method is the Time Dependent Hartree-Fock (TDHF) equation for the 
one-body density matrix p T (r 1 ,r 2 ,t) = (r 1 |p T (t)|r 2 ) : 



■ h dp T 



where H T is the one-body self-consistent mean field Hamiltonian depending implicitly on 
the density matrix and r is an isotopic spin index. It is convenient to modify equation 
(1) introducing the Wigner transform of the density matrix 



/s s 
d 3 s exp(-ip • s/h)p T (r + - , r - - , t) 

and of the Hamiltonian 



(2) 



H T w (r 



, p) = J d 3 s exp(— ip ■ s/h)(r + 



r--)- 



(3) 



Using (2,3) one arrives [3] at 



or 



■ sin 



(v)Mvy - (v^ • (v/H i^r, 



(4) 



dt h { 2 

where the upper index on the bracket stands for the function on which the operator in 
these brackets acts. It is shown in [5, 6], that by integrating equation (4) over the phase 
space {p, r} with the weights x^x^ . . . Xi k pi k+1 . . ■Pi n _ 1 Pi n , where k runs from to n, one 
can obtain a closed finite set of dynamical equations for Cartesian tensors of the rank 
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n. Taking linear combinations of these equations one is able to represent them through 
irreducible tensors, which play the role of collective variables of the problem. However, 
it is more convenient to derive the dynamical equations directly for irreducible tensors 
using the technique of tensor products [7]. For this it is necessary to rewrite the Wigner 
function equation (4) in terms of cyclic variables 

(>J! | sin I? E (-1)° [(V-X • (VS) / - {Via)* ■ (Va)'l I Hfrr, (5) 



dt h 2 „ . 

1 a=—l 



with 



1,0 . d . d I d . d 

1 r ■ ^ 1 / • x 

and the analogous definitions for V+ x , Vq , , andp + i , p , The required 

equations are obtained by integrating (5) with different tensor products of r a and p a . Here 
we consider the case n — 2. 



2.1 Model Hamiltonian, Equations of motion 

The microscopic Hamiltonian of the model is 

A -2 1 2 z v 

^ — 1 — 2 £ j 

y 2 Z N 

+ 2 K E (-1HE (fc-^rO^fo) + E ^(r^fo)}, (6) 

where the quadrupole operator g 2/i = ^16^/5 r 2 F 2 ^ and iV, zJ are the numbers of neutrons 
and protons, respectively. The mean field potential for protons (or neutrons) is 

V T (r,t) = ImwV + E (-l)^2%W^(r), (7) 

where = kQ^ + , = + and the quadrupole moments Ql^(t) 

are defined as 

Ql li (t) = Jd{p,r}q 2lx (r)r(r,p,t) 

with fd{p, r} = 2(2irh)~ 3 jd 3 p jd 3 r, where the factor 2 appears due to summation over 
spin degrees of freedom. To simplify notation we omit spin indices, because we consider 
spin saturated system without the spin-orbit interaction. 
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Substituting spherical functions by tensor products r 2 Y 2fM = \ —r\ u , where 

I Off ^ 

and Cxa % \ v is the Clebsch-Gordan coefficient, one has 

V^\mu^^- i r z l-A- (8) 

Here 

z^^xR^ + xR^, z^ = x R^ + xR n 2^ x = 6«, x = 6«, 

^(0 = /4P,rK/(r )P ,t). (9) 

Integration of equation (5) with the weights r\ , {j'p) x ^ = {r <S> p}\^ and p\ yields 
the following set of equations [1]: 



-W -- 

it v m 

^ - ^ + - 2>/5 £ ^i + lffl} Wa*. = 0, A = 0,1, 2 



-rR\u L xa = °, A = 0,2 

dt v m X/x 



^ + 2m^-4v^Ev / ^{2l{}(^I)A, = 0, A = 0,2 (10) 

where {^l} is the Wigner 6j-symbol and the following notation is introduced 

Px,(t) = Jd{p,r}plj T (v,p,t), L^(t) = Jd{p,v}(rp) x ,r(r,p,t). 

By definition q 2)l = x/Qr^, Q T 2fi = V6R^, R^ = -Ql /V3 with Q T m = N T < r 2 > being 
the mean square radius of neutrons or protons. The tensor L\ v is connected with angular 
momentum by the relations L[ = ^7^3, L\ ±1 = |(7J T 

We rewrite equations (10) in terms of the isoscalar and isovector variables = R x ^ + 
R v Xfl1 R XfJi = R x — R\ (and so on) with the isoscalar kq = (k + R)/2 and isovector K\ = 
(k — R)/2 strength constants. There is no problem to solve these equations numerically. 
However, we want to simplify the situation as much as possible to get the results in 
analytical form giving us a maximum of insight into the nature of the modes. 

1) We consider the problem in small-amplitude approximation. Writing all variables 
as a sum of their equilibrium value plus a small deviation 

Rx,(t) = RH + n Xll (t), P x ,(t) = PH + V^(t), L x ,(t)=LH + C x ,(t), 



we linearize the equations of motion in 1Z\^, and 7Z\^, P\^, C\^. 

2) We study non-rotating nuclei, i.e. nuclei with LH = L{1 = 0. 

3) Only axially symmetric nuclei with R 2 ±2 = R%±i = -^2±2 = -^2±i = are considered. 

4) Finally we take 

^8 = ^ = 0. (n) 

This means that equilibrium deformation and mean square radius of neutrons are supposed 
to be equal to that of protons. 

Due to the approximation (11) the equations for isoscalar and isovector systems are 
decoupled. Further, due to the axial symmetry the angular momentum projection is a 
good quantum number. As a result, every set of equations splits into five independent 
subsets with ji = 0, ±1, ±2. The detailed derivation of formulae for eigenfrequencies and 
transition probabilities together with all necessary explanations are given in [1]. Here 
we write out only the final results required for the comparison with respective results 
obtained in the framework of RPA. 

2.2 Isoscalar eigenfrequencies 

The isoscalar subset of equations with /j, — 1 is 
K 21 - 2C 21 /m = 0, 

C21 - V 2 i/m + [muo 2 + 2« (<22o + 2Qoo)] ^21 = 0, 
V21 + 2[muj 2 + K Q e 2 l]jC 21 = 0, 

£11 = 0. (12) 

Imposing the time evolution via e~ tnt for all variables one transforms (12) into a set of 
algebraic equations. The eigenfrequencies are found from its characteristic equation which 



reads 



- 4c 2 - ^(QH + \qH)\ = 0. (13) 
m 3 

moj 2 _ 2 u 2 

For hq we take the self-consistent value = — — — — , where u = (see Appendix A) 

4Qoo 1 + 3d 

4 

with the standard definition of the deformation parameter Q20 — Qoott Then 

o 

Q 2 [Q 2 -2u 2 (1 + 5/3)] = 0. (14) 



The nontrivial solution of this equation gives the frequency of the n — 1 branch of the 
isoscalar GQR 

Q 2 = Ql = 2u> 2 {\ + 5/3). (15) 

Taking into account the relation (A. 7) we find that this result coincides with that of [9]. 
The trivial solution Q = Q = is characteristic of nonvibrational mode corresponding 
to the obvious integral of motion C u = const responsible for the rotational degree of 
freedom. This is usually called the 'spurious' or 'Goldstone' mode. 

2.3 Isovector eigenfrequencies 

The information about the scissors mode is contained in the subset of isovector equations 
with /i = 1 

£21 - 2C 2 i/m = 0, 

C 2 i-V2i/m+ m u 2 + kQII + 4kiQII TZ 2 i — 0, 
V21 + 2[mu 2 + - 6K0Q20 Ai = 0, 

£11 + 3RQ e 2 lK 2 i = 0. (16) 

Imposing the time evolution via e~ %nt one transforms (16) into a set of algebraic equations. 
Again the eigenfrequencies are found from the characteristic equation which reads 



9 ^fi 

fi 4 - n 2 [Aoo 2 + - Kl Q% + — («i + 2k )Q2] + ^(«o - «i)«o(<52 9 o) 2 = 0. (17) 
m m m z 



Supposing, as usual, the isovector constant K\ to be proportional to the isoscalar one, 
Ki = an , and taking the self-consistent value for k , we finally obtain 

ft 4 - 2Q 2 u 2 (2 -a)(l + 5/3) + 4c^ 4 (l - a)5 2 = 0. (18) 

The solutions of this equation are 



n 2 ± = Gj\2 - a)(l + 5/3) ± sju>\2 - a) 2 (l + 5/3) 2 - 4D 4 (1 - a)5 2 . (19) 

The high-lying solution f2 + gives the frequency Q iv of the fi — 1 branch of the isovector 
GQR. The low-lying solution f2_ gives the frequency Q sc of the scissors mode. 
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We adjust a from the fact that the IVGQR is experimentally known to lie practically 
at twice the energy of the isoscalar GQR. In our model the experimental situation is 
satisfied by a = — 2. Then 



^ = 4^(l + ^ + ^/(l + ^-^ 2 j, tf„ = tf\\ + *--^ + ^-\p ) - (20) 

2.4 Linear response and transition probabilities 

A direct way of calculating the reduced transition probabilities is provided by the theory 
of the linear response of a system to a weak external field 

Fit) = Fexp(-iftt) + F t exp(iftt), 

where F = J2f=i fs is a one-body operator. A convenient form of the response theory is 
e.g. given by Lane [8] (see also section 4). The matrix elements of the operator F obey 
the relation 

| < u \p\0 > | 2 = h lim (ft - Q u )< ifj\F\ip > exp(-iftt), (21) 

where |0 > and \v > are the stationary wave functions of the ground and unperturbed 

excited states; ip is the perturbed wavefunction of the ground state, VL V = (E v — E )/h are 

the normal frequencies, the bar means averaging over a time interval much larger than 

1/ft, ft being the frequency of the external field F(t). 

Magnetic excitations 

z * eh 
F = K = E A/M A/. = -^V(ry lM ) • [r x V]^ = 7 (rp)i, p N = — . (22) 



2mc V 2tt 

Ft(MU -91 <- nr\F*>\()> \2_^~^ m ^ n ,2 ^c-2(l + V3)^ 2 2 

u(juf \ \ -01 ^ ? - ? ,.|FP|n> I 2 - ^"^n 2 ^^^-2(1 + ^/3)^ 2 

These two formulae can be joined into one expression by the simple transformation of the 
denominators. Really, we have from (19) 

±(ft 2 ,-ft 2 c ) = ±(ft 2 + - ft 2 ) = ±2^(2 - a) 2 (l + 6/3y - 4^(1 - a )p 

= 2ft|-2cj 2 (2-a)(l + 5/3) = 2ft|-(2-a)(cj 2 + w 2 ). (25) 
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Using these relations in formulae (23) and (24), we obtain the expression for the B(M1) 
values valid for both excitations 



, l-amu 2 2 n 2 u - 2(1 + 5/3)u 2 



B(M1)„ = 2| < HFfilO > I' = — — O^ Vw-^-akl + f/S)] ^ <26) 
Electric excitations 



F = H> = T,fM> k = er 2 Y 2 , = (3rl, P = eJ^. (27) 
< ^\F^ >= PR% = \(5{R 2 , - R 2 »). 



B(F2\ 21 r,,ipp|n. l 2 - 6 ^ 5 r> (i + V3)^ e -2(^) 2 , . 

5(£2) sc - 2| < SC |F 21 |0 > | - — -Qoo sUSn^m • ( } 



Rf/^ -2l< ? d^ln>l 2 - e2/l5 o (i + V3)^-2(^) 2 

5(£2) ro _ 2| < ,t,|F 21 |0 > | - — -Qoo nUnl _ nl) • (29) 



B{E2) is = 2| < i a |j&|0 > | 2 = ^^Qoo[(l + */3)n£, - 2(^) 2 ]/[fiJ 3 . (30) 

Using relations (25) in formulae (28) and (29) we obtain the expression for the B(E2) 
values valid for all three excitations 

B(F2) -2\<v\Fno>\*- e2h 5 O a + V3)^-2(^) 2 

B(E2) V -2\< u\F 21 \0 > | - — Y^Qoo nAn2 _- 2{2 _ a){1 + s/3)y (31) 

The isoscalar value (30) is obtained by assuming a — 1. 



3 Random Phase Approximation (RPA) 

Standard RPA equations in the notation of [3] are 

^ ] {^ijmn-^-nj [^ij^mn i^m Q) ~l~ ^inmj] Yfij} ^-^^rni- (32) 
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According to the definition of the schematic model by [3], the matrix elements of the 
residual interaction corresponding to the Hamiltonian (6) are 

f, .. _ K , n-r* T)t' 

with D im =< i\q2i\m > and /t nn = k pp = k, K np = R. This interaction distinguishes 
between protons and neutrons, so we have to introduce the isospin indices r, r' into the 
set of RPA equations (32): 

( e m - e i) X mi + K rT'Dl^D T jn Xlj + ^ K TT > DJ^D^jY^ = HflX^, 

n,j,r' n,j,r' 

E "tt>DZD&X* + {e m - e[]Y^ + £ ^D^D^ = -/Kll^. (33) 

n,j,r' n,j,r' 

The solution is 



J~)T* D T * 

x -' = w^, KT ' Y ^ = -mf^r (34) 



With e mi = e m - C I aIld RT = £r' ^rr'C T 

The constant C T is defined as C T = Y,n,j(Dj n X^ j + D^ lj Y^ j ). Using here the expressions 
for X^j and Y^ given above, one derives the useful relation 

C T = 2S T K T = 2S T «Vr'C T ', (35) 

t' 

where the following notation is introduced: 

S T = Y,\DU 2 F2 e }L )2 (36) 

mi V t mi/ 

with E = HQ. Let us write out the relation (35) in detail 

C n - 2S n («;C n + RC P ) = 0, 

C p - 2S p (RC n + kC p ) = 0. (37) 

The condition for existence of a nontrivial solution of this set of equations gives the secular 
equation 

(1-2S ii k)(1-2S p k) -AS n S p R 2 = 0. (38) 

Making obvious linear combinations of the two equations in (37), we write them in terms 
of isoscalar and isovector constants C = C n + C p , C = C n — C p 

C - 2(S n + S p )k C - 2(S n - S p )k x C = 0, 

C - 2{S n - S p )k q C - 2(S n + S p )k x C = 0. (39) 
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Approximation (11) allows us to decouple the equations for isoscalar and isovector con- 
stants. Really, in this case S n = S p = 5/2; hence, we obtain two secular equations 

1 - 2Sk = 0, or 1 - Sk = Sk (40) 

in the isoscalar case and 

1 - 2Sk! = 0, or 1-Sk= -Sk (41) 

in the isovector one, the difference between them being in the strength constants only. 
Having in mind the relation K\ = an Q , we come to the conclusion that it is sufficient to 
analyze the isovector case only - the results for isoscalar one are obtained by assuming 
a — 1. 



3.1 Eigenfrequencies 

The detailed expression for the isovector secular equation is 



AK 1 mi ^ e r, 



-mi 



The operator D has only two types of nonzero matrix elements D mi in the deformed 
oscillator basis. Matrix elements of the first type couple states of the same major shell. 
All corresponding transition energies are degenerate: e m — €{ = h(u x — oo z ) = eo- Matrix 
elements of the second type couple states of the different major shells with AiV = 2. 
All corresponding transition energies are degenerate too: e rn — e« = h{uj x + uj z ) = e 2 . 
Therefore, the secular equation can be rewritten as 



I e D : e 2 D 2 

2 



2«i E^-d + E^-e 2 2 (43) 



The sums D — ^ \D mi \ 2 and D 2 — ^ \D m i\ 2 can be calculated analytically 

mi(AN=0) mi(AN=2) 

(see Appendix B): 

D = — ? e , D 2 = —^e 2 . (44) 
Let us transform the secular equation (43) in polynomial form 

E A - E 2 [(e 2 + 4) + 2 Kl (e D + e 2 D 2 )} + [e 2 e 2 2 + 2^^^ + e 2 D )} = 0. 
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Using here the expressions (44) for D , D 2 and the self-consistent value of the strength 
constant (A. 3), we find 



E 4 - E\l - a/2)(e 2 + e\) + (1 - a)e 2 e 2 = 0, 

or 

n 4 - n 2 (2 - a)u 2 + + (1 - a)ut = 0, (45) 

with the notation u\ = u 2 + oj 2 and cut = iyo 2 — oj 2 ) 2 . This result coincides with that of 
[2]. By a trivial rearrangement of the terms in (45) one obtains the useful relation 

n 2 (n 2 -uj 2 + ) = (i- a )(n 2 uj 2 + -uji). (46) 

Inserting expressions (A. 3) for oj 2 , uj 2 into (45), we find uj\ = 2lj 2 (1 + 5/3), u 4 = 45 2 u 4 
and reproduce formula (18) for the isovector case 

VI 4 - 2fiV(2 - a)(l + 5/3) + 4co 4 (l - a)5 2 = 0. 

Taking here a = 1 we reproduce formula (14) for the isoscalar case 

tt 4 -2tt 2 tu 2 (l + 5/3) = 0. 



3.2 B(E2)-factors 

A 

According to [3], the transition probability for the one-body operator F = ^ f s is Calcu- 
li 

lated by means of the formulae 

< o|f> >= Y,(fLx% + fX), < H^ r |o >= E(/™^r + (47) 

mi mi 

Quadrupole excitations are described by the operator (27) with f 2fl = er 2 Y 2fl = eD, where 
e = t\j y|^. The expressions for X^, Y mi are given by formulae (34). Combining these 
results we get 

< 0|F 2 > >= TeKl £ \DU \ 2 '7p r 2 = 2eff*5* = eC». (48) 



mi 



The constant is determined by the normalization condition 

°v,v' — V mi mi 1 mi 1 mi )i 
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that gives 



(Cu) 2 



I n p \ 2 f p (r n v i n 



mi I 



{siy [ei - ( e L) 2 ] 2 (cv p ) 2 (s»y [El - (e^ t ) 2 } 2 

The ratio C n /C p is determined by any of the equations (37): 

C n 1 - 2S p k 2S D R 



(49) 



(50) 



Cp 2SPR 1 - 2S*k 
Formula (49) is considerably simplified by the approximation (11), when S p = S n = 
5/2, e^j = ej^j, -D]^ = D^. Applying the second forms of formulae (40,41) it is easy to 
find that in this case C n /C p = ±1. As a result, the final expression for B(E2) value is 

B{E2) V = 2| < 0|Fj> > | 2 = 2e~ 2 f 16i^ 2 £ \D mi \ 2 " m \ ) . (51) 

With the help of formulae (44) this expression can be transformed into 

B(E2) V 



5 e 2 Q 



oo 



8tt muj 2 a 2 E v 



e 2 e 2 

fc i fc 2 



(E 2 u -e 2 ) 2 (E 2 u -e 2 ^ 2 



5 e 2 Q 00 (£ 2 -e 2 ) 2 (E 2 -e 2 ) 2 
8tt mu 2 a 2 E u {El - e|) 2 eg + (E 2 - e 2 ) 2 e 2 2 
5 e 2 /iQoo {n 2 u uo 2 + -uji) 2 



(52) 



16tt mu 2 !], _ + W 2 W 4 ■ 

At first sight, this expression has nothing in common with (31). Nevertheless, it can be 
shown that they are identical. To this end, we analyze carefully the denominator of the 
last expression in (52). Summing it with the secular equation (45) (multiplied by u\ ), 
which obviously does not change its value, we find after elementary combinations 

Denom = fi^ 2 - 2Q 2 u ujt + w 2 c/ + u 2 + [ft* - ft 2 (2 - a)uj 2 + + (1 - a)u 4 _] 
= iv 2 + n 2 u [2n 2 u - (2 - a)u 2 + ] - ui[2Ql - (2 - a)u 2 + ] 
= (nlcol -c/)[2fi 2 -(2-«V 2 ]. (53) 

This result allows us to write the final expression 



i 



B(E2) V = -^-— 2 Qoo Q 21 . ( 54 ) 

16% mu 2 \l u [2\V u — (2 — a)uj+\ 

which coincides with (31) (we recall that u 2 ^ = 2u 2 (l + 5/3), uj a _ = 45 2 u 4 ). By the simple 

transformations this formula is reduced to the result of Hamamoto and Nazarewicz [2] 

5 e 2 h 

(taking into account, that they published it without the constant factor Qoo)- 

62nmu 
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3.3 B(Ml)-factors 



In accordance with formulae (22), (47), (34) the magnetic transition matrix element is 

(fll)imDf m (fll)miD 



<o\Fr 1 \u>=K^Y: 



)P* 

mi 



(55) 



E - f p . E + e p ■ 
As it is shown in Appendix B, the matrix element (/n)j m is proportional to D p m (formula 
(B.16). So, expression (55) is reduced to 



< OlFftli/ > 



eft 



= A? 



eft 



D p D p * 

^ im im 



D p D p *- 

^ mi mi 



p 

rn i ) 



u cV5 



(^-^) P ^E^ 



e p -LE 2 

L m! L i/ 



P \2]' 



(56) 



With the help of approximation (11) and the expressions (44) for D , D 2 we find 



lP i , CI eh . 2 2 Qoo / 

: (^-^J77— o( 



<0|F 1 p 1 |z/> = 



2^ Cv /5 



:y5 



mw 2 a^f^u; 2 , — u;i) 



(57) 



2 cv^o" 21 zy ^ 

Relation (46) and the self-consistent value of the strength constant K\ = ano were used 
in the last step. For the magnetic transition probability we have 

u 4 (1-a) 2 



B{Ml\ 



2| < Ol-F&lz/ > | 2 = 2 



{CD 2 e~ 2 , A (1 - a) 



-£(£2),. (58) 



4 5c 2 - ft 2 20c 2 ill 
This relation between £?(M1) and B(E2) was also found (up to the factor l/(20c 2 )) 
by Hamamoto and Nazarewicz [2]. Substituting expression (54) for B(E2) into (58) we 
reproduce (with the help of relation (46)) formula (26). 



3.4 "Synthetic" scissors and spurious state 

The nature of collective excitations calculated with the method of Wigner function mo- 
ments is quite easily revealed analyzing the roles of collective variables describing the 
phenomenon. The solution of this problem in the RPA approach is not so obvious. That 
is why the nature of the low-lying states has often been established by considering overlaps 
of these states with the "pure scissors state" [11, 12] or "synthetic state" [2] produced by 
the action of the scissors operator 

s x = n-\< if > i p x - < if >4 n ) 
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on the ground state 

\Syn >= S x \0 > . 

In the considered model the overlap of the "synthetic" state with the real scissors mode 
(and with IVGQR) can be calculated analytically. Surprisingly, it was not done until now. 
Let us at first modify the definition of the "synthetic" state. Due to axial symmetry one 
can use the I* component instead of I£, or any of their linear combinations, for example, 
the \i = 1 component of the magnetic operator F{ , which is much more convenient for us. 
The terms < I t2 > are introduced to ensure the orthogonality of the synthetic scissors to 
the spurious state \Sp >= (I n +I p )\0 >. However, we do not need these terms because the 
collective states \v > of our model are already orthogonal to \Sp > (see below); hence, the 
overlaps < Syn\u > will be free from any admixtures of \Sp >. So, we use the following 
definitions of the synthetic and spurious states: 

\Syn >= AT 1 ^ - Fl\)\0 >, \Sp >= (Ffi + F^)\0 > . 

Let us demonstrate the orthogonality of the spurious state to all the rest of the states 
\v >. As the first step it is necessary to show that the secular equation (38) has the 
solution E — 0. We need the expression for S T (E = 0) = S T (0). In accordance with (36), 
we have 

eoA) 



S T (E) 



+ 



e 2 D 2 



£2_ £ 2 E 2_ e 2 



S T (0) = - 



CO ^2 



The expressions for Dq, D t 2 are easily extracted from formulae (B.10), (B.ll): 



-^o ~~ *3oo 
m 



1 + p 



D 2 — Qoo 
m 



1 + 15 



So we find 



(59) 



S T {0) 



-—Q T 
m 



1 + 15,1 



+ -) + 

£2 eo 



1-15,1 



e 2 



eo' 



h 2 A5 T Ql 



oo 



3g 



20 



=2 e 



Wl(w2- W 2)T' 



where, in accordance with (B.12), 
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m 



(60) 



(61) 



Finally, we get 



2S P (0) 



Q p 



20 



20 



15 



2^(0) = -—%—, 1-2S*(0)«- 



20 



It is easy to see that substituting these expressions into (38) we obtain an identity; there- 
fore, the secular equation has a zero energy solution. 

For the second step it is necessary to calculate the overlap < Sp\v >. Summing (56) 
with an analogous expression for neutrons, we get 

L V " r mi mi\ v c miJ 

iH ,E V £ KM - "IY E 77, \ D }£% w - (62) 



cv/5 u ^ uy x ZJ (4i) T (i£ " eL) 1 
Applying the algebraical identity 

1 J_A _L 



e 2 (£ 2 -e 2 ) £ 2 V £ 2 -e 2 ' 
and remembering the definition (36) of S T we rewrite (62) as 

eh 



eh Kl 



In accordance with (35) and (50), 



ui - coins- - sp(o)) + (coi - coins- - s«m§ 



.(63) 



^ = l ~ 2SVK (64) 

Noting now (see formula (60)) that (cj 2 — co 2 z ) T S T (Q) = —^Ql and taking into account 
relations (61), we find 

( 1—2 9 p k 

<Sp\u> = (3\[(kQI + kQI)2S v -Ql] + [{KQl + RQl)2S n -Ql\ 



2S n R 

i _ 2 9 p k 

^ [(2S P « - l)Ql + 2S P /^] + [(2S n K - + 2ff n /cQg)] 2Sn _ 

1-2S p k 



2S*RQl + (25V - 1)Q5 2ffl _ | 



= {2,S n /t25 p /t - (1 - 25 ,n /c)(l -2,S p /t)} = 0, (65) 

3 eh K p 

where (3 = 7= — - and Q2 = Q20 The expression in the last curly brackets coincides 

m cv 5 £i, 

obviously with the secular equation (38) that proves the orthogonality of the spurious state 
to all physical states of the considered model. So we can conclude that strictly speaking 
this is not a spurious state, but one of the exact eigenstates of the model corresponding to 
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the integral of motion I n + I p . In other words [3]: "In fact these excitations are not really 
spurious, but they represent a different type of motion which has to be treated separately." 
The same conclusion was made by N. Lo Iudice [13] who solved this problem approximately 
with the help of several assumptions (a small deformation limit, for example). 

The problem of the "spurious" state being solved, the calculation of the overlaps 
< Syn\v > becomes trivial. Really, we have shown that < OlF^ + F^u >= 0. That 
means that < 0\F^\u >= - < 0|Fft|i/ >; hence, < Syn\u >= A/" 1 < 0|F& - F^\u >= 
< OlFftli/ > and 

U 2 = \< Syn\v > | 2 = 2^- 2 B(Ml) u . (66) 

The nontrivial part of the problem is the calculation of the normalization factor A/". It is 
important not to forget about the time dependence of the synthetic state which should 
be determined by the external field: 

\Syn(t) >= AT'KFfi - F^)e^ nt + (F[\ - F^e^O > . 

Then we have 

N 2 = 2<0|(F 1 p 1 -F 1 n i) t (A P i-A n i)|0> 

= 2£ < 0\(F?i ~ F^Ylph >< phKPr, - F^)\0 >= 2]T | < p/iKFfi - F^)\0 > \ 2 

ph ph 

= 2Y / \<ph\F[ 1 \0>\ 2 = 2Y / \(fn)ph\ 2 - (67) 

r,ph T,ph 

With the help of relation (B.16) we find 

<ph\r 2 Y 21 \0 > | 2 V 



r,ph \ 



t ph 



Expressions for Dq, D\, uj t x , uo t z are given by formulae (59), (B.12). To get a definite 
number, it is necessary to make some assumption concerning the relation between neutron 
and proton equilibrium characteristics. As usual, we apply the approximation (11), i.e., 
suppose Qoo — Qoo) Q20 = Qm- ^ ls eas y to check that in this case formulae for are 
reduced to the ones for the isoscalar case, namely (A. 3), and Dq = D /2, D T 2 = D 2 /2, 
where D and D 2 are given by (44). So we get 



87T 2c mu 2 Ve e 2 J 2n k 
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The estimation of the overlap for 156 Gd with 5 = 0.27 gives Af 2 = 34.72/4 and U 2 = 
0.53 (see eq. (66)), that is two times larger than the result of [11] obtained in QRPA 
calculations with the Skyrme forces. The disagreement can naturally be attributed to the 
difference in forces and especially to the lack of pair correlations in our approach. In a 
small deformation limit U 2 = \\J\ ~ 0.6. 

This is the maximum possible overlap of the "pure" (or "synthetic") scissors with 
the real scissors. The increasing of 5 and /or taking into account pairing correlations 
decreases its value, that is confirmed by numerous microscopic calculations with various 
forces [17]. Such small overlap leads inevitably to the conclusion, that the original model 
of counter rotating rigid rotors [10] has not very much in common with the real scissors 
mode, the correct description of which requires the proper treatment of the Fermi surface 
deformation and the coupling with IVGQR. 

3.4.1 Superdeformation 

A certain drawback of our approach is that, so far, we have not included the superfluidity 
into our description. Nevertheless, our formulae (20, 26) can be successfully used for the 
description of superdeformed nuclei where the pairing is very weak [2, 10]. For example, 
applying them to the superdeformed nucleus 152 Dy (5 ~ 0.6, huo = 41/A 1//3 MeV), we get 

E iv = 20.8 MeV, B(Ml) iv = 15.9 /4 

for the isovector GQR and 

£ sc = 4.7MeV, B(M1) SC = 20.0/4 

for the scissors mode. There are not so many results of other calculations to compare 
with. As a matter of fact, there are only two papers considering this problem. 
The phenomenological TRM model [10] predicts 

E iv ~26MeV, B{Ml) iv ~ 26/4, £ sc ^6.1MeV, B{M1) SC ~ 22 fi 2 N . 

The only existing microscopic calculation [2] in the framework of QRPA with separable 
forces gives 

S iw ~28MeV, B(Ml) iv ~ 37/4, £ sc ^5-6MeV, B(M1) 1+ ~ 23 /4- 

18 



Here B(M1) 1 + denotes the total Ml orbital strength carried by the calculated K w = 1 + 
QRPA excitations modes in the energy region below 20 MeV. 

It is easy to see that in the case of IVGQR one can speak, at least, about qualita- 
tive agreement. Our results for E sc and B(M1) SC are in good agreement with that of 
phenomenological model and with E sc and B(M1) 1 + of Hamamoto and Nazarewicz. 

It is possible to extract from the histogram of [2] the value of the overlap of calculated 
low-lying 1 + excitations with the synthetic scissors state: | < Syn\l + > | 2 w 0.4. The 
result of our calculation U 2 = 0.43 agrees with it very well. So the natural conclusion of 
this section is that the correct treatment of pair correlations is obligatory for a reasonable 
description of the scissors mode. 

3.5 Equations of motion 

Let us look on WFM equations of motion from the RPA point of view. Is it possible to 
construct something similar in the RPA approach? Equations (10) are written for average 
values of operators and are valid for the description of the arbitrary amplitude motion. 
One can compare with RPA only their linearised version, obtained by the variation of 
equations. The variables of linearised equations are the variations of the above mentioned 
average values. It is natural to suppose some correspondence between the variation of 
the average value of F operator and the matrix element of the type < 0\F\u > used 
to calculate transition probabilities. To check this idea we have to derive dynamical 
equations for matrix elements of the operators r\ , p\ and {rp)\^ to compare them with 
linearised equations (10). To this end we combine RPA equations (33) in accordance with 
the definition (47) of matrix elements: 

hQ, v y](fim-^mi + fmiXmi ) = e mi(fi m ^-m,i ~ fmiXmi ) + ^ ~ fmiXmi)- (^) 

mi mi mi 

Taking into account the relations 

^mifim [ft -^o]im> ^mifmi [ft -^o]mi> 

one rewrites this equation as 

an„ < o|f> >= J2{[f\ Hi\ im x^ + \f, + K{fj m DZ - f mi DZ)}- (7i) 

mi 



19 



The Hamiltonian of the axially deformed harmonic oscillator corresponding to the mean 
field (8) is 

#o r (r) = £{f^ + - 2 m^vl + Zl {eq)rUs)}. (72) 

Let us consider the operator / = \/§r 2l = <?2i = D. Calculating the commutator 

2 

[r 21 ,H ] = ih—(rp) 21 
m 

we find from (71) the following equation 

«i„<o|5»> = tftVe- E{((rp)2i)Sn« + ((rp)2i)M"} 

s=l m mi 

J2(DJ m DZ - D^DZ). (73) 

mi 

Taking into account relations (D*) im = (D)* mi and |-D m j| 2 = |-Dj m | 2 we find, that the last 

sum in (73) is equal to zero. Applying again formula (47) and introducing the notation 

A ^ A 

R\n = XX r sW> L ^ = E( r «P«W we write (73) as 

s=l s=l 

2 

< 0\R T 21 \v >= — < 0\L T 21 \v > . (74) 

Identifying the matrix elements < Ol-R^ik > an d < 1 Z/^ | ^ > with 1Z 2l and £^1 respec- 
tively we reproduce the variation of the first equation in (10) (having in mind the time 
dependence via e~ lQt ). 

Let us consider the operator / = {rp) 2 i. The required commutator is evaluated to be 

h h 
[(rp) 2 i, H ] = i—p 2 2l - iTimuj 2 r 2 2l - i-=Z 20 (eq)rl v 

With this result equation (71) looks as 

hQ u < 0\L 2l \v >= i— < 0\P 21 \u > -ihmcu 2 < OlR^u > 

-i\zi,{eq) < 0|i&|i/ > +KlY\{{rp) T 21 )imDZ - ((r^U^Cl, (75) 
V6 mi 



where the notation P Xfl = XX^sW has been introduced. The last sum is calculated with 



s=l 

the help of formula (B.17). Using the fact, that e im = —e mi , one gets 

E[(M 2 i)L^ - ((rp) 2i y mi DZ] = -i^ E tmi[(r 2 2i y im DZ + (r 2 2i y mt DZ] 

mi mi 

rji rji A 

= -ij-^T,z T mi\D T mi \ 2 = "^(W + ^D 2 ) = -^^(1 + 5/3)Ql . 
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According to the definitions (see formula 35) we have 



Kl = 5>n T C£ = ( X < 0\R n 2l \u >+x< 0\R p 21 \u >)/VQ, 

T 

K = Y.^rCl = (x < 0\R p 2l \u >+ X < 0\R n 21 \u >)/v / 6. 

T 

So, the equation (75) (let us say, for neutrons) is transformed into 
-i£l v < 0\L 21 \u >= +^ < Ol-RJili/ > -mu 2 < 0\R 2l \u > 

-±=Z» Q {eq) < 0\R^\u > - 2 -{l + 5/3)QUx < 0\R^\u > +X < 0\R p 21 \u >). (76) 

The equation for protons is obtained by interchanging indices n and p. One has to compare 
this equation with the variation of the second equation in (10) with A = 2, fj, — 1. Let us 
write this variation in detail: 

|^ 1 -i-n i +ma; 2 ^ 1 -2v^ £ yfij + lffi} £ C% Jv [Zl,(eq)K%+6ZLir jv (eq)] = 0. 
ai m j=o,2 

We recall, that only Rl (eq) and R 2Q (eq) have non zero values, so this equation is reduced 
to 



-C T 21 - -V T 21 +mcu 2 n T 21 - mi¥i}cil 21 z^(eq)n T 21 



d „ r 1 

at m 

-2v / 55^ 1 [{^}C 2 2 1 1 i00 ^ (eg) + >fi$&}C%~K^eq)\ = 0. 



In agreement with definition (8) of Z\ its variation is 

sz^ = x n% + xn\ , sz^ = x n% + xK, . 

Substituting 6j-symbols and Clebsch-Gordan coefficients by their numerical values we 
obtain finally (e.g. for neutrons) 

4*^ " -Ki + mu 2 -R} 2 \ + ±Z» (eq)K» 21 + + 5/?>)QUxKi + X^i) = 0. 
at m y/Q 6 

This equation coincides obviously with (76) if to assume the time dependence via e~ tnt 
and to identify the matrix elements < Q\R 2l \v >, < OlL^ik > an d < 1 ^21 1 ^ > with the 
variables 1Z 21 , C 21 and V 21 , respectively. 

Let us consider the operator / = (j3 2 ) 2 i. The required commutator is 

2 

l(p 2 ) 21 , H ] = -ih2muj 2 (rp) 21 + iMy/E^ p3 + HlHyC^Z^eq^rp)^ 

21 



and one obtains from (71) the following equation 

h£l v < Q\p2\\v > = -ihlmuj 2 < Q\L T 21 \ v > 

+ HlHyClljAieq) < 0|L> > 

+K;J2{((P 2 )2ihmDZ - miiUDZ]- (77) 

mi 

It is easy to show (with the help of formula (B.18) that the last sum is equal to zero. 
This equation must be compared with the variation of the last equation in (10) with 
A = 2, /i = 1. Let us write it in detail. Taking into account that L T XfM (eq) = we find the 
equation 

+ 1mu?C7 21 - 4v/5E pi + l^CH^Zl^C], = 
at j=1 

that obviously coincides with (77) if to assume the e~ tnt time dependence and to identify 
the proper RPA matrix elements with the respective WFM variables. 

4 WFM versus RPA 

The exact relation between RPA matrix elements and the respective WFM variables can 
be established with the help of the linear responce theory. Let us first recall, following 
Appendix D of [3], the necessary definitions concerning the density and the density matrix. 
The density operator is defined as 

A 

P( r ) = 5^( r - r s) =Y, d Pi( r ) a l a ^ (7 8 ) 

s=l pq 

where d pq (r) =< p\S(r — r)\q >= ^ < P*p{ rcrT )4>q{ v<TT ) an d 4>q{ raT ) are single-particle wave 

<7,T 

functions. Indices p, q include spin and isospin quantum numbers a and r. 

The density of particles in the system depends on its state \I/ and is defined as the 
average value of a density operator over this state: 

p(r) =< tt|p(r)|W >=Y,d Pq (r)p qp 

pq 

= A J2 I d 3 r 2 .-d 3 r A \^(rar, r 2 cr 2 r 2 , r A a A r A )\ 2 , (79) 

a,T,...,a A ,TA 
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where p qp =< ^\a) p a q YS> > . The particle density (79) can be interpreted as the diagonal 
element (in the coordinate space representation) of the density matrix which is defined 

as 

p(rar,r'a'r') = 5>;(rW)0,(r<7r) < *\a},a q \* >= £ d M (rW, var)p qp (80) 

pq pq 

with d pq (r'a'r' ,rar) = <p*(r'a'r')(f) q (rar). The average value of the arbitrary one-body 
operator 

F = Y.fs = Y.UA a * ( 81 ) 

s=l P<? 

is written in terms of the density matrix as 

< >= fpq < *l<»Kl* >= E /m/°9P = Tr (/P)- 

P<3 P9 

Let us consider the system to be in the weak external time-dependent field 

W(t) = Wexp(-int) + W*exp(int), (82) 

where W = J2 pq w pq a l a q is a one-body operator. The change of the ground state wave 
function produced by this field is found by using the time- dependent perturbation theory 
[14]: 

^(t) = \0> +J2W> [cve~ tnt - c* u e lUt ] . (83) 

V 

Here |0 > and \v > are stationary eigenstates of the unperturbed system and 

< u\W\0 > _ „ < u\a p a q \0 > _ < 0\W\u > _ „ < Ola+ajz/ > 

(84) 

Inserting this expression into formula (80) we obtain the perturbed density matrix 

p{rar, r'a'r' , t) = po(rar, r'a'r') + 5p{rar, rVV, t), 
where poirar, r'a'r') is the unperturbed (equilibrium) density matrix 

p (rar, r'a'r') = ^ ^(rW, rar) < 0|aja p |0 >= ^ ^(rW, rtrr)^ 

P<? P<3 

and Sp(rar, r'a'r', t) is the change of the density matrix 

5p(rar, r'a'r', t) = £ <V(rVr' , r<7r)p« (t) (85) 

P<3 
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with 



int 



Ppqtt) = E [(< 0|aja p |^ > c u - < z/|aja p |0 > c v )e 

+ (< v\a\a p \® > c*- < 0|aja p |z/ > c*)e im ] . (86) 

Deriving (85) we neglected the terms proportional to \W\ 2 . At this stage it is necessary 
to remind that we work in a Hartree-Fock approximation. That means that stationary 
states |0 >, \v > are Slater determinants; matrix p p °J = p q S pq is diagonal with p q — 1 for 
levels below the Fermi level and p q = for levels above the Fermi level. The requirement 
(p + dp) 2 = (Po + bp) leads to the well known [3] property of the matrix p^: it has only 
particle-hole nonvanishing matrix elements. Looking to formula (86) we see that it is 
possible for the matrix elements < 0|aja p |z/ > to be different from zero only for particle- 
hole combinations of indices q,p. Consequently, the summation over p, q in formula (84) 
for c n and c n will also be restricted only to particle-hole pairs. So we can write p p )J as 



P$(t) = E [R Pq ,M^ nt + Rl P M^) e 

p'q' 

where 



W p i q i , 



Rpq,p>q>(ty - E 



< 0|aja p |z/ >< z/|ap,a g /|0 > < 0|ap,a ? /|f >< z/|a^a p |0 > 



is the RPA response function [3], where the index pairs pq and p'q' are restricted to 
particle-hole pairs. For the change of the arbitrary operator average value we have: 

5<^\F\^>=J2f Pq P qP ) - (87) 

pq 

Now we are ready to analyze the WFM variables. The first one is 



R^(t) = 2(2nh)- 3 J d 3 p Jd 3 rrlr(r,p,t). 



Using here the definition (2) of the Wigner function and the definition of the 5-function 
we find 

Rl^t) = j^^Jd'rr^Jd'sJd'pexpi-tp-s/^p^r + ^r-^t) 
= 2 [d 3 r r^p T (r, r, f ) = £ f d 3 r r 2 Xlx p(rar, rar, t) 
= EE [(Prrl^KraT^iraT) < V\a],a q \V > 



pq a 
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pq pq 

= <*IE^( a )l*>=<*l*W*>- 

s=l 

For the variation of this variable one can write the following chain of relations 
SRl,(t) = ^(t) = 2|rf 3 rr^(r,r,t) 



= EE d*rrl^{TOT)U™T)pU{t) = 

pq cr pg 

= E(< 0|^a> > c v - < v\Rl^ > c v )e~ m 

V 

+ E(< > < I^> > 3)***- 

V 

For the second variable we have 

L^(t) = 2(2^ Jd 3 p Jd 3 r(rp) x ,r(r,p,t) 

= (2nh) 3 I d?>r I I ^ exp(_ip ' s / h )P T ( r + |> r - |> *) 

= ^{(rlV-VVlry,*)}^ 
= -yE /^{MV-V'D^KrVr.t)}^ 
-yEE /rf 3 r{0;(rar)(rV) AM g (r^) 

Z pg <r J 



(? (r ( 7r)(rV) AM 0;(r ( 7r)} < ^|aja 9 |^ > 

/S 
2 



E{(Mam)p T , + tt^fco E / ^ 0;(rar)0 9 ( rC rr)}p gp (t) 



= < *| EMs)amI^ > 5aoA^ t =< > +ih— 5 X0 N T . 



s=l 

The variation of this variable is 



SL\,(t) = ^(t) = -ft/^{(r[V-V']) A ^(r,r',f)} f 

= E((^)am)p>Sp } W 
py 

v 



The third variable is 



P^(t) = 2(2^)- 3 |rf 3 p Jd 3 rplj T (r,p,t) 
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= (2^)3 / ^ / rf3fi / ^ pp %> ex PHP ' s /% T ( r + |, r - |, 
= -| /rf 3 r{(V-V')^(r,r',t)} r=r , 

= "T E E W""0 V^^(r^r) + ^(ror) V^;(r<rr) 

-2[V0 g (nxr)V0;(nxr)] AM } < ^aja^ > 
= "^EE l<Pr(t>l{raT)Vl^ q {r<JT) < V\a)p, q \V > 

= ECpW^W =< *l EpK s )I* >=< MH^ > • (92) 

pq s=l 

The variation of this variable is 

= ^(^-/^{(V-V'^r/.t)}^ 

pg 

= E[<0|A>>c,-<H^|0>c,]e-^ 

V 

+ E(< HAVIO > <£- < 0|P A > > *C]e** (93) 

V 

The structure of variables Hxn, ^ AAt , Pa^ (89,91,93) demonstrates in an obvious way the 
relation between the WFM method and RPA. One sees, for example, that the dynamical 
equations for the WFM variables Ttx^ is a linear combination of the dynamical equations 
(74) for the transition matrix elements < OI-Ra/J^ >, the mixing coefficients c n and c n 
being determined by the structure of the wave packet (83). Naturally, the same is true 
for the variables Cx^, Vxp,- The dynamical equation for < 0|i?A^|^ > is in turn, the linear 
combination of RPA equations (33) (or (32) in the case of arbitrary interaction) for the 
amplitudes X pq , Y pq , the mixing coefficients being particle-hole matrix elements of the 
operator Rx^. 

As we see, there exist exact relations between the dynamical equations for the variables 
of the WFM method (moments) and the RPA dynamical equations for the amplitudes 
X pq , Y pq . One should note however, that these relations are exact only in our simplified 
model, because in general both methods, to be exact, have to operate within an infinite 
number of dynamical equations. In RPA one replaces the infinite number of particle-hole 
pairs of the shell model by the infinite number of phonons with the hope that the essential 
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part of physics is described by the small number of the lowest energy collective phonons 
and, consequently, one can neglect the rest of (infinite number) phonons. The coupling 
of the dynamical equations for X pq , Y pq , corresponding to different particle-hole pairs is 
realized by the matrix elements of the nucleon-nucleon interaction (see equations (32)). 
An analogous situation is observed in the WFM method, where the dynamical equations 
for Cartesian tensors of rank n = 2 are coupled (by the interaction terms in (5)) with 
dynamical equations for tensors of rank n — 3, these equations being coupled with the ones 
for tensors of rank n = 4 and so on up to n = oo. And again one hopes that the essential 
part of physics is described by a few number of the lowest ranks tensors. This hope is 
based on the evident consideration that the higher rank tensors (moments) are responsible 
for more refined detailes and that, by neglecting them, one does not appreciably influence 
the description of the more global physics which is described with the lower ranks tensors. 
In this game of including only the lowest rank tensors one has to remember the trivial 
(but important) rule: the highest rank of tensors must not be less than the multipolarity 
of the studied motion. 

It is easy to see that the nature of truncation in the two methods is quite different. So 
that in practical calculations with realistic Hamiltonians one can not establish the exact 
relation between these methods unless one works in the full space in both methods. 

Of course there are exceptions like the case of the mean field potentials with quadratic 
coordinate dependence (harmonic oscillator with quadrupole-quadrupole or monopole- 
monopole residual interaction). Due to the huge degeneracy of the particle-hole configu- 
ration space all RPA sums are calculated analytically without any approximations. The 
same happens in WFM method - the dynamical equations for tensors of different ranks 
decouple and one obtains a finite set of equations, which can be solved exactly. As a con- 
sequence, both methods give identical results for integral characteristics of the collective 
motion, such as energies and transition probabilities. 

A difference appears in the description of various distributions in coordinate space, for 
example, transition densities and currents, where the WFM method can not give the exact 
result, because it deals only with integrals over the whole phase space {p,r}. However, 
in principle the WFM method can give any number of moments required to construct 
approximate expressions for these distributions (see below). 
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4.1 Flows 



We want to know the trajectories of infinitesimal displacements of neutrons and protons 
during their vibrational motion (the lines of currents). The infinitesimal displacements 
are determined by the magnitudes and directions of the nucleon velocities u(r, t), given 
by 

/Acfip 
^-^p/(r, p,t) 

4 



^3 / rf3pPeXp (~ ip ' S /^^( r+ ^' r ~^' t ) 



ifi 



-2tk{(V - V)p(r, r', t)} r=r , = -- E{(V - V)p{var, r'ar, t)} r= 

ih 

-yEE{^HV^ 9 (H - <f) q (rar)V(j);(rar)} < ^|aja g |^ > 

pq tr,T 

m Y.h q ij)p q p{'t) = m < | E 3 P q( r Wp a <i\^ > 



pq 



pq 



m < *|J(r)|* > . 



(94) 



The current density operator j(r) has the standard quantum mechanical definition [3]: 



A t A 

Ar) = E^W = 7T- £[*(r - f.)V. + V s 5(r - r s )] = Y.UM^ 

s=l Z7TO s=l P9 

ft 



j OT (r) =< p| — [<f(r - f)V + V5(r - r)]|g > 
^ ^[^ rar ) V< ^ r(Tr ) " ^(«"W;(«7t)] = 4^[0;(r)V0 g (r) - ^(r)V^(r)]. 



2mi 



The variation of u generated by the external field (82) is 
p e «(v)6u(v,t) = £j P9 (r)p«(t) 

pq 

= £[< 0\J(t)\v > c u - < u\J(t)\0 > c v \e~ m 

v 

+ £[< H^(r)|0 > cl- < 0\J(v)\u > cj] e «*. 



(95) 



To proceed further two ways are possible. 

The first, so to say direct way is obvious. Having solutions (34) for X^, Y^ d we can 
calculate transition currents with the help of formula (47): 

f D* i D* ■ 1 



< 0\J(v)\u >= J2(j im X^+j mi Yr) = KJ2 

mi mi 

"?' D* i D* 

Jtrn^tm Jmi^mi 



E 



v mi(A7V=0) 



E u — e E v + e 



+ E 

mi(AN=2) 



Ely ^mi E v -\- 6j 
-E v — €2 E u + €2 



(96) 
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The operator D has a finite number of particle-hole matrix elements D mi , so, in principle, 

the sums in (96) can be calculated exactly. The same is true for the coefficients c v . 

Therefore, one can find the exact (in the frame of RPA) result for the velocity distribution 

5u(r, t). However, even in this simple model one can not find a compact analytical 

expression for sums in (96) - the field of velocities can be constructed only numerically. 

The second way allows one to derive an approximate analytical expression for Su(r, t). 

The main idea lies in the parametrization of the infinitesimal displacements £ (r, t) [5] . 

dtAr, t) 

Let us recall the main points. By definition <5iij(r,t) = — — . The displacement is 

at 

parametrized [1] by the expansion 

3 3 3 

&(r,i) =G t (t) +J2 G iA t ) x j+ X) G i,jk(t)xjX k + G i,jki(t)xjX k xi H (97) 

j=l j,k=l j,k,l=l 

which, in principle, is infinite, however one makes the approximation keeping only the first 
terms and neglecting all the rest of it. For example, in [1] only the two first terms were 
kept. It turned out, that Gi = due to the triplanar symmetry of considered nuclei. The 
coefficients Gij were expressed analytically in terms of the variables 72-21 (t) and Cu(t). 
Using the dynamical relations between 72-21 (t) and Cu(t) given by the last equation of 
the set (16), the final formulae for £j(r, t) were found to be 

6 = V2BJ 13 x 3 , 6 = V2BJ 23 x 3 , £ 3 = V2A(J 13Xl + J 23 x 2 ) 

with 

Jia = (72 2 -i - 72 2 i)/2, J23 = KK2-1 + 72 2 i)/2, 
A = ^[l-2^{l-a)5]/[Q m {l- 2 -5)l 

B = ^[1 + 2^(1 - a)5]/[Q 00 (l + p)}, (98) 

The pole structure of the right hand side of equation (95) tells us, that the transition 
current can be calculated by means of an expression analogous to (21): 



< 0\Ji(r)\u >= h lim (fi - 0,)p eg (r)&(r, t) exp(^t)/ < u\W\0 > . (99) 

— SZi/ 

For the & from above we obtain (using formulae (89) and (84)) 

B 

< 0|Ji(r)|i/ >= -iVL u p eq {r)—= < 0\R 2 -i - R 2 i\v > x 3 , 

V2 
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< 0\J 2 (r)\u >= n u p eq (r)—= < 0\R 2 ^ + R 2l \v > x 3 , 

v2 

A - 
< 0| J 3 (r)\u >= -iVL v p eq (r)—=[< 0\R 2 -i - R 2 i\v> x l + i< 0|-R 2 -i + #21 |f > x 2 \. 

v2 

It is obvious that the second way is more adequate for the WFM method, because the 
moments 1Z 2 \(t) and C n (t) are just WFM variables and the dynamical relation between 
them is just given by the WFM dynamical equation. 

If necessary, one can find the next term of the series (97). To calculate the respective 
coefficients Gijki{t) in the WFM method one is obliged to derive (and solve) the set of 
dynamical equations for higher (fourth) order moments of the Wigner function. To solve 
the same problem with RPA, it is necessary to construct the analogous set of dynamical 
equations for transition matrix elements of the respective operators. The required work 
is approximately of the same order of complexity in both cases. 

In conclusion in full RPA one must calculate the currents numerically leading to fine 
detailes (shell effects) whereas in WFM and approximate RPA treatment one obtains their 
gross structure with analytical formulas. The latter feature is quite important in order to 
understand the real character of the motion under study since current patterns produced 
numerically from complicated formulas with a lot of summations like in (96) can hardly 
be interpreted physically. A good example is the interplay of the scissors mode and the 
isovector giant quadrupole resonance. Looking only at the flow patterns (see Figs. 1, 2 in 
[1]) one would not be able to tell that the former is mostly rotational with a small amount 
of an irrotational component and the other way round for the latter, as this can be seen 
from eqs. (42)-(47) in [1]. In this respect it is important to work with the infinitesimal 
displacements £j, because by definition they are differentials (£i = dx,!; 2 = dy,^ = dz) 
which allow one to construct differential equations for the current fields. For example 
equation (99), showing that transition current is proportional to a differential, allows 
one to derive a differential equation for the current field in RPA directly from (96). For 
example 




<0|J 2 (r)|z/> dy 
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5 Conclusion 



The properties of collective excitations (the scissors mode, isovector and isoscalar gi- 
ant quadrupole resonances) of the harmonic oscillator Hamiltonian with the quadrupole- 
quadrupole residual interaction (HO+ QQ) have been studied with two methods: WFM 
and RPA. We have found that both methods give the same analytical expressions for 
energies and transition probabilities of all considered excitations. This, however, does not 
mean that WFM and RPA are identical approaches in all respects. For example current 
distributions are described differently in the two approaches even in this simple model. In 
general both methods are not equivalent unless one makes sure that the space of moments 
corresponds exactly to the particle-hole space used in RPA. However, the spirit of WMF 
is rather to drastically reduce the dimensions in considering only low order rank tensors. 
In this way, of course, one will loose the fine structure in the spectrum but still the gross 
structure will be well approximated. One also may check the convergence of the method 
in increasing the number of moments. In the case of well defined resonances only some 
more satellites to the main peak should appear. Such a method may be particularly useful 
in the case of deformed nuclei where the dimension of the RPA matrices becomes easily 
prohibitive. 

It makes no sense to speak about advantages or disadvantages of one of the two 
discussed methods - they are complementary. Of course, RPA gives complete, exhaustive 
information concerning the microscopic (particle-hole) structure of collective excitations. 
However, sometimes considerable additional effort is required to understand their physical 
nature. On the contrary, WFM method gives direct information on the physical nature 
of the excitations. Our results serve as a very good illustration of this situation. What 
do we learn about the scissors mode and IVGQR from each of the two methods? RPA 
says that the scissors mode is mostly created by AiV = particle-hole excitations with a 
small admixture of AiV = 2 particle-hole excitations and vice versa for IVGQR. Without 
further effort - this is about all. One does even not suspect the key role of the relative 
angular momentum in the creation of the scissors mode. On the other hand, the WFM 
method directly reveals that the scissors mode appears due to oscillations of the relative 
angular momentum with a small admixture of the quadrupole mode and vice versa for 
IVGQR. Further, it informs us about the extremely important role of the Fermi surface 
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deformation in the formation of the scissors mode. This demonstrates very well the 
difference between two approaches: the RPA describes the fine structure of collective 
excitations whereas the WFM method gives their gross structure. 

Two new mathematical results are obtained for the HO+QQ model. We have proved 
exactly, without any approximations, the orthogonality of the "spurious" state to all phys- 
ical states. In this sense, we have generalized the result of Lo Iudice [13] derived in a small 
deformation approximation. The analytical expressions are derived for the normalization 
factor of the synthetic scissors state and overlaps of this state with eigenstates of the 
model. It is shown, that the overlap of the synthetic scissors with the real scissors reaches 
its maximal value ~ 0.6 in a small deformation limit. The increasing of 5 and /or taking 
into account pairing correlations decreases the overlap, that is confirmed by numerous 
microscopic calculations with various forces [17]. Such small overlap leads inevitably to 
the conclusion, that the original model of counter rotating rigid rotors [10] has not very 
much in common with the real scissors mode, the correct description of which requires 
the proper treatment of the Fermi surface deformation and the coupling with IVGQR. 

Appendix A 

It is known that the deformed harmonic oscillator Hamiltonian can be obtained in a 
Hartree approximation "by making the assumption that the isoscalar part of the QQ 
force builds the one-body container well" [15]. In our case it is obtained quite easily by 
summing the expressions for V p and V n (formula (7)): 

V(v,t) = \(V*(T,t) + V n (r,t)) = ImcuV + ^o £ (-l)"Q 2 -„(*)?2„(r). (A.1) 

In the state of equilibrium (i.e., in the absence of an external field) Q 2 ±i = Q2±2 = 0. 
Using the definition [16] Q 2 o = Qoof ^ and the formula g 20 = 2z 2 — x 2 — y 2 we obtain the 
potential of the anisotropic harmonic oscillator 

nr) = yk(* 2 + rWV] 

with oscillator frequencies 

u 2 x = uj 2 y =uj 2 {l + a5), u 2 = u 2 {l-2a5), 
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where a = —k q — — The definition of the deformation parameter 5 must be reproduced 
3mur 

by the harmonic oscillator wave functions, which allows one to fix the value of a. We have 

Qoo = -— + — + —, Q 20 = 2-— -— , 

m oo x uj y uo z m uo z uo x 

where E x = T<f =1 (n x + -)j and n x is the oscillator quantum number. Using the self- 
consistency condition [16] 

where S and u are defined in the spherical case, we get 

Q20 _ uj 2 x - u 2 z 2aS 4 
Q 00 ul + 2ul I -ad 3 ' 

Solving the last equation with respect to a, we find 

< A ' 2 » 

Therefore, the oscillator frequences and the strength constant can be written as 

col = co 2 y = co 2 (l + ^), co 2 z =co 2 (l- 2 -5), ko = ~ (A.3) 

with Co 2 = u) 2 /{l + The condition for volume conservation uj x uj y uj z = const = lOq 
makes 00 (^-dependent 

co 2 -co 2 1 + 5* 

CO — LOn 



(1 + |^) 2 / 3 (1 - f^Vs- 
So the final expressions for oscillator frequences are 

^^(i±l) ,/3 , ^(^|) 2/3 . (A . 4) 

It is interesting to compare these expressions with the very popular [16, 3] parametriza- 
tion 

co 2 x = co 2 = co' 2 (l + 2 -5'), u 2 z =co' 2 (l-p'). 
The volume conservation condition gives 



12 _ w o 



CO = 



;i + |5') 2 / 3 (i- 15') 1/3 ' 

so the final expressions for oscillator frequences are 



,2 _ , ,2 _ , ,2 x 3 1 , ,2 _ , ,2 I 3J 



«:=^=^(ir^j • w --= w HttfJ ■ (A - 5) 
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The direct comparison of expressions (A. 4) and (A. 5) allows one to esatablish the following 
relation between 5 and 5': 



1 + 25' 1-25'' 
One more parametrization of oscillator frequences can be found in the review [17]: 



, ,»2 , ,5)2 



2 2 ^ 2 
LV„ = UJ„ = t— , UJ 



x y l-f<5"' 2 1 + 
One has from the volume conservation condition 

^" 2 = c 2 (l-^") 2/3 (l + ^") 1/3 , 
so the final expressions for oscillator frequences are 

that coincide exactly with (A. 4), i.e. <5" = 5. 

It is easy to see that equations (A. 4) correspond to the case when the deformed density 
p(r) is obtained from the spherical density po( r ) by the scale transformation [9] 

(x,y,z) — > (xe a/2 ,ye a/2 , ze~ a ) 

with 

which conserves the volume and does not destroy the self-consistency, because the density 
and potential are transformed in the same way. 

It is necessary to note that Qoo also depends on 5 

V00 = — ( 1 1 ) = — SoUM - 2 ^ 2> = ^OO/i , 4ni/vi 2 r\2/V 

where Qq = A^R 2 , R = r Q A l / z . As a result, the final expression for the strength constant 
becomes 

mm 2 ( I- fdV /3 _ mu 2 

4Q8o U + I^y 4g§ ' 

that coincides with the respective result of [9]. 
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Appendix B 



To calculate the sums D = ^2 \D m i\ 2 and D 2 = ^2 \D m i\ 2 we employ the sum- 

mi(AN=0) mi(AN=2) 

rule techniques of Suzuki and Rowe [9] . The well known harmonic oscillator relations 



n x +l) i 



Px^n x = ~i \ m ^ x ( y/n^i/j nx _ i - V^ + li.+i) (B.l) 



allow us to write 



XZl\) nx l\) nz = — = ( y/n x n z tl) nx - l^„, - 1 + J (n x + l)(n z + 1)^+1^+1 

+ yj (n x + l)n 2 ^+i^n z -i + + IjVV-iVv+i), 

PxPz ^ / 

/ * VW0n* = ~o }= ( \f n xn z i>n x - 1 - 1 + J (n x + l){n g + l)lp nx+1 1p nz+1 

m z u x uj z 2m^u x uj z v 



-yj (jl x + - V n ^( n ^ + 1 )^n.,-l^n,+l)- ( B - 2 ) 

These formulae demonstrate in an obvious way that the operators 

P — \{ zx H r~ and P 2 = ^(zx p^) 

contribute only to the excitation of the AiV = and AiV = 2 states, respectively. Fol- 
lowing [9], we express the zx component of r 2 Y 2 \ = \J j^D = — ^J^z(x + iy) as 

zx = P + P 2 . 

Hence, we have 

A A 

e o \ < ^\J2 z sX s \mi > \ 2 = e J2\ < 0\^2P (s)\mi > \ 2 

mi(AN=0) 8=1 mi s=l 

1 <0|EP (*),[ff,i;Po( a )]]|0>, (B.3) 



2 ,, i ! 



where eo = h(u x — uj z ). The above commutator is easily evaluated for the Hamiltonian 
with the potential (A.l), as 

< 0|EP ( S ), [h£p o ( 8 )]]\0 >= -^o f <0|Ell ^ |0> - <0|Ell ^ |0> ) • (B.4) 

8=1 s=l Zm \ V x Ld z ) 
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Taking into account the axial symmetry and using the definitions 

a ^4 
Qoo =< 0| E(2^ + ^ 2 )|0 >, Q20 = 2 < 0| ^(4 ~ ^)|0 >, Q20 = Qoo-S, 

s=l s=l 6 

we transform this expression to 

A A h /1 + -5 l--5\ 

< 0|E Po(s), [H, E P o (s)]]\0 >= -%Q o 8- . (B.5) 

s=l s=l 6m V U z J 

With the help of the self-consistent expressions for u> x , uj z (A. 3) one comes to the following 
result: 

< oignw, p.gfl.WP >= |f| = (J - 1) 2 . (b,) 

By using the fact that the matrix elements for the zy component of r 2 Y" 2 i are identical to 
those for the zx component, because of axial symmetry, we finally obtain 

<0L r /2lF> = T7. ~ € = Ta 2 i h) ■ B ' 7 

mi(AAr=o) *=i 16nmu 2 16vr m u 2 \l - ^ J 

By calculating a double commutator for the P 2 operator, we find 

2 ™S rf , 1 0| S r ' l2l|,ra > 1 - T^^' 2 " T^^=? ll^f*J ' (&8) 

where e 2 = h(cu x + uj z ). 

We need also the sums Dq and D 2 calculated separately for neutron and proton systems 
with the mean fields V n and V p , respectively The necessary formulae are easily derivable 
from the already obtained results. There are no any reasons to require the fulfillment 
of the self-consistency conditions for neutrons and protons separately, so one has to use 
formula (B.5). The trivial change of notation gives 

z z h fl + -5 p l--5 p \ 

< OlEPoOO, [# p ,£Po( S )]]|0 >= f-egggo £— , (B.9) 

8 = 1 8 = 1 bm \ UJX UJz J 

Z <=> h /l _|_ ±XP 1 _ 2 rp\ 

4 E I<oie^ 2 iK>I 2 = t|--«o [-^--^- , (B.10) 
% E I<oiE^iK>I 2 = t^-^o i± |^ + ^^ • (b.ii) 

™ ( ^= 2 ) s =i 167rm V ^ ^ / 
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The nontrivial information is contained in oscillator frequences of the mean fields V p and 
V n (formula (7)) 



«) 2 = c 2 [i - -^(«g3o + «gs,)], (c 2 n ) 2 = c 2 [i + -^(«c^o + (B.12) 



The above-written formulae can be used also to calculate the analogous sums for 
operators containing various combinations of momenta and coordinates, for example, 
components of an angular momentum, tensor products (rp) 2 i and (p 2 ) 2 i- Really, by 
definition 1\ = yp z — zp y , I 2 = zp x — xp z . In accordance with (B.l), we have 

h 



xp z ^n x ^n z = ~i 77 J — ( \Jn x n z i\) nx -xi\) nz _ 1 - J(n x + l)(n z + l)ip nx+1 ip. 
z v uo 



n z +l 



+yj (n x + lKW+iVv.-i - ^n x (n g + (B.13) 



Therefore, 



^«in z = - \ / — ) ( \/n x n z il) nx -iip nz -i ~ \l{n x + l){n z + 



2 V V w 



cu 2 



— + \l—)(\f( n x + l ) n z^n x +i^n z ~i - Jn x (n z + (B.14) 
z y y 

Having formulae (B.2) and (B.14), one derives the following expressions for matrix ele- 
ments coupling the ground state with AiV = 2 and AiV = excitations: 



< n x + 1, n z + 1\I 2 \0 >= i 



.fi{u? x -u? z ) (n x + l)(n z + l) 



2 oo x + 00 z 



< n x + l,n z — 1\I 2 \0 >= i 



(n x + l)n z 



2 L0 X — 00 z 



00 X U0 Z 



, 11 m & K + l)(n 2 + 1) 
2m V 



< n x + l,n z — l\xz\0 >- 



h (n x + l)n z 



2m V 



It is easy to see that 



(to 2 - to 2 ) 

<n x + l,n z + 1 \I 2 >= im^^ — < n x + l,n z + lbz O >, 

+ ^z 

(uo 2 — cj 2 ) 

< n x + l,n z - 1|/ 2 |0 >= im— — < n x + l,n z — l\xz\0 > . 

OJr - CO* 



(B.15) 
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Due to the degeneracy of the model all particle-hole excitations with AN = 2 have the 
same energy e 2 and all particle-hole excitations with AN = have the energy e . This 
fact allows one to join the last two formulae into one general expression 

(u 2 — u 2 ) 

< ph\I 2 \0 >= ihm— — < ph\xz\0 > . 

Taking into account the axial symmetry we can write the analogous formula for I±: 

(ijj 2 — lu 2 ) 

< ph\li\0 >= -ihm— — < ph\yz\0 > . 

The magnetic transition operator (22) is proportional to the angular momentum: /i±i = 

ie I 3 - 

\ — (I 2 =F Hi) Therefore, we can write 

4mc V 27r 

< pfc|/i±i|0 >= eM^-"* 2 ) K ph \ r i Y2±l \0 > . (B.16) 

2cV5 e P A 

Similar calculations for the tensor product (rp) 2 i = —\[zp x + xp z + i(zp y + yp z )} lead to 
the following relation: 

Tfl / 2"7r TYl 

<ph\(rp) 2 i\0 >=i — \ —e ph < ph\r 2 Y 2±1 \0 >= i^e^ < ph\r 2 2l \0 > . (B.17) 
n \ id in 

Two kinds of particle-hole matrix elements are obtained from the second formula of 
(B.2): 

. . \{n x + l)(n z + 1) 

<n x + l,n z + l\p x pz\0 >= -hmuj x uj z \ — , 

V 2u x 2u z 

1|A A , / (n x + l)n z 

<n x + l,n z - l\p x p z \0 >= hmuj x uj z \ — — . 

V 2u x 2u z 

Simple comparison with (B.15) shows that 

< n x + 1, n z + l\p x p z \0 >= —m 2 uj x uj z <n x + l,n z + l\xz\0 >, 

< n x + 1, n z — l\p x p z \0 >= m 2 uj x uj z < n x + 1, n z — l\xz\0 > . 
With the help of obvious relations 

2u x u z = uj 2 x + lj 2 - e 2 /h 2 , -2uj x u z = uj 2 £ + uj 2 - e 2 2 /h 2 
these two formulae can be joined into one expression 

o 

Tfl 

< ph\p x p z \0 >= + <*>l - e 2 ph /h 2 ) < ph\xz\0 > . 
By definition p| x = —p z {p x + ip y ) and r\ x = —z[x + iy), hence, 

2 

Tfl 

< Ph\p 2 21 \0 >= —(co 2 x + ^1 - e 2 ph /h 2 ) < ph\T 2 21 \0 > . (B.18) 



38 



References 

[1] E. B. Balbutsev, P. Schuck, Nucl. Phys. A720 (2003) 293; 

E. B. Balbutsev, P. Schuck, Nucl. Phys. A728 (2003) 471, Erratum 

[2] I. Hamamoto, W. Nazarewicz, Phys. Lett. B297 (1992) 25 

[3] P. Ring and P. Schuck, The Nuclear Many Body Problem (Springer, Berlin, 1980) 

[4] S. Chandrasekhar, Ellipsoidal Figures of Equilibrium (Yale University Press, New 
Haven, Conn. 1969) 

[5] E. B. Balbutsev, Sov. J. Part. Nucl. 22 (1991) 159 

[6] E. B. Balbutsev, P. Schuck, Nucl. Phys. A652 (1999) 221 

[7] D. A. Varshalovitch, A. N. Moskalev and V. K. Khersonski, Quantum Theory of 
Angular Momentum (Nauka, Leningrad, 1975) 

[8] A. M. Lane, Nuclear Theory, (Benjamin, New York, 1964) 

[9] T. Suzuki, D. J. Rowe, Nucl. Phys. A289 (1977) 461 

[10] N. Lo Iudice, La Rivista del Nuovo Cimento 23 (2000) N.9 

[11] R. R. Hilton et al., Proc. 1st Int. Spring Seminar on Nuclear Physics "Microscopic 
Approaches to Nuclear Structure" (Sorrento, 1986) ed A. Covello (Bologna: Physical 
Society), p. 357 

[12] A. E. L. Dieperink, E. Moya de Guerra, Phys. Lett. B189 (1987) 267 
[13] N. Lo Iudice, Nucl. Phys. A605 (1996) 61 

[14] L. D. Landau and E. M. Lifshits, Quantum Mechanics (Nauka, Moscow, 1974) 
[15] R. R. Hilton, Ann. Phys. (NY.) 214 (1992) 258 

[16] A. Bohr and B. Mottelson, Nuclear Structure, vol. 2 (Benjamin, New York, 1975) 
[17] D. Zawischa, J. Phys. G: Nucl. Part. Phys. 24 (1998) 683 



39 



